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ABSTRACT 

There are very few confirmed black holes with a mass that could be ~ 4 M 0 and no neutron stars 
with masses greater than ~ 2 Mq , creating a gap in the observed distribution of compact star masses. 

Some black holes with masses between 2 and 4 Mq might be hiding among other X-ray sources, whose 
masses are difficult to measure. We present new high-speed optical photometry of the low-mass X- 
ray binary V1408 Aql (= 4U 1957+115), which is a persistent X-ray source thought to contain a 
black hole. The optical light curve of V1408 Aql shows a nearly sinusoidal modulation at the orbital 
period of the system superimposed on large night-to-night variations in mean intensity. We combined 
the new photometry with previously-published photometry to derive a more precise orbital period, 

P = 0.388893(3) d, and to better define the orbital light curve and night-to-night variations. The 
orbital light curve agrees well with a model in which the modulation is caused entirely by the changing 
aspect of the heated face of the secondary star. The lack of eclipses rules out orbital inclinations greater 
than 65°. Our best models for the orbital light curve favor inclinations near 13° and black hole masses 
near 3M 0 with a 90% upper bound of 6.2 Mq, and a lower bound of 2.0 Mq imposed solely by the 
maximum mass of neutron stars. We favor a black hole primary over a neutron star primary based on 
evidence from the X-ray spectra, the high spin of the compact object, and the fact that a type I X-ray 
burst has not been observed for this system. Although uncertainties in the data and the models allow 
higher masses, possibly much higher masses, the compact star in V1408 Aql is a viable candidate for 
a black hole lying in the mass gap. 

Subject headings: binary systems: X-Ray binaries, individual (4U 1957+115) 


1. INTRODUCTION 

Low mass X-ray binaries (LMXB) are comprised of a 
compact stellar remnant, either a black hole or a neu¬ 
tron star, that accretes matter from a low-mass Roche- 
Lobe-filling secondary star via a circumstellar disk. The 
LMXB 4U 1957+115 was first detected by the Uhuru 
satellite (Giacconi et al. 1974) and its o ptical counter¬ 
part V 1408 Aql was first identified by Margon et al. 
(1978]). The X-ray properties of 4U 1957+115 are un¬ 
usual. LMXBs usually cycle between active and quies¬ 
cent states, but 4U 1957+115 has been persistently active 
for more than 40 years, the longest interval of any known 
LMXB, remaining in a spectrally soft, disk-domin ated 
X-ray state with no detect able radio jet (Wijnands et al. 
2002 Russell et al. 2011). V1408 Aquiiae may be the 
only BH LMXB, with the possible exception of LMC X- 
3, that has been found persistently active. We note that 
LMC X-l, LMC X-3, and Cyg X-l ( |Russell et aT~|[2CiTo| ) 
are persistent high mass X-ray binary systems. 

The X-ray light cu rve of the system has n ot shown any 
orbital modulat ion ([Wijnands et a h1|2002|), but optical 
observations by Thorstensen "[ 1987) revealed a nearly si¬ 
nusoidal orbital variation with a peak-to-peak ampli tude 
of 23% and a period of 9.33 hr. Hakala et al. (1999) ob¬ 
served the light curve on two nights and saw a change 
in the shape of the light curve and a significant increase 
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in the amplitude of the variations. Further evidence of 
night-to-night variations was presented by|R.ussell et al.| 


(20101 from 144 images of the source spread over three 


years. The larger data sets of Bayless et al. (2011) and 


Mason et al. (2012) confirmed both the sinusoidal orbital 


modulation and the chan ges in mean brightn e ss from 


night to n ight. Following Thorstensen (1987), Bayless 


et al. (20111 showed that the orbital light curve can be 


reproduced by a model in which the secondary star is 
heated by flux from the accretion disk. The orbital mod¬ 
ulation is produced entirely by the heated face of the 
secondary star as it rotates into and out of view . Our 


results differ from the models by Hakala et al. (2014) 
in that their models do not find any evidence for a sec¬ 
ondary star in the spectral energy distribution, and do 
not agree with a secondary whose face is irradiated by 
X-rays. The SED measured in this work is consistent 
with our model. Since the heated face of the secondary 
would be just a small, high-temperature black-body per¬ 
turbation on the SED and would not be distinguishable 
in the measured SED. 

Several lines of evidence indicate that 4U 1957+115 has 
a black hole primary. It is common for LXMBs with neu¬ 
tron star primaries to show type I X-ray bursts but 4U 
1957+115 has never shown such bursts. The X-ray spec¬ 
trum is well described by a multi-temperature blackbody, 
with an additional non-thermal power law com ponent in 


15% of the observations (Nowak et al. 2012). Its high 
inner disk temperature and small in ner disk radius are 


1983; Wijnands et al. 2002 

Nowak et a T. 2008; Kussell et| 

al. 

2010 

Nowak et ai.||2012 

). Models of the X-ray spec- 
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tra that allow for rotation of the primary star yield large 
spin rates: from a* > 0.9 for a black hole with a mass of 
3 Mq and a distance of fO Kpc, to a near-maximal spin of 
a* « 1 for larger values of the mass and distance (Nowak 


et al. 2012). The observed spins of black holes range up 


to a’ ~ 1, while the spi ns of neutron star s are exclu¬ 
sive ly less than a* = 0.1 (Mi ller et al. ||2011 ). Yaqoob et 


all dl993|), |Singh et al.| (|l994|) , and lllayless et al.| Q2011J 
argued that the primary in 4U 1957+115 is a neutron 
star, in part because the mass ratio q = M 2 /M\ is large, 
suggesting a small-mass primary that is more consistent 
with a neutron star than a black hole. We note, however, 
that a black hole with an unusually low mass could also 
yield a high mass ratio. 

The distribution of the known black hole and neu¬ 
tron star masses has a gap between 2 Mq and ~ 4 Mq. 
Among the neutron stars with reliably measured masses, 
the most massive is J0348+0432 at 2.01 ± 0.21, M ,m fol- 


lowed by J1614+2230 at 1. 97 ± 0.04 (Demo rest et 

ij |Latimer| 

0422+32 h 


al. 2010 Antoniadis et al. 2013 
least massive black hole, (IRC) J 


3.97 ± 0.95 Mq (Gelino & Harrison 


2UT4+ The 

has a mass of 
2003), followed by 


GRS 1009-45 with a mass of near 4)4 +1+ and not less 
than 3.64M 0 (Filippenko et al. 19991, and poss ibly by 


4U 1547-4 7 with an estimated mass of < 4 Mq (Kreid- 


berg 20121. A very massive neutron star or a very low 


mass black hole could be pr oduced fr om a progenitor 
with a mass of ~ 22 Mq (Fryer 1999). The probabil¬ 
ity that the observed gap between the masses of neutron 


stars and black holes is a n 

lere statistical fluke is low 

(Ozel et al. 

2010 

Farr et al. 

2011 

1, although it remains 


possible that the gap has a non-zero but sparse popula¬ 


tion. 


If the mass gap is real, it has important implic ations 


for th e physics of core-co l lapse supernovae (Fryer et al. 


2012|. Belczynski et al. (20121 proposed a theoretical 


explanation tor the existence of the mass gap that de¬ 
pends on the growth time of instabilities that lead to 
core-collapse supernovae. Stars in the 20 — 40 Mq are 
the ones likely to have high mass neutron stars or low 
mass black holes as remnants that would lie in the mass 
gap. If the growth time of instabilities is larger than 
200 milliseconds, a continuous distribution would be ex¬ 
pected in this range. But models that assume a growth 
time of 10 - 20 mill iseconds introduce and explain the ob¬ 
served mass gap (Belczynski et al. 2012). The observed 
mass distribution might, however, be subject to strong 
selection effects. One possibility is that low-mass black 
holes are hiding among other X-ray sources whose masses 
are notoriously difficult to measure. The absorption-line 
spectra o f t heir se condary stars are generally not visible. 


Ozel et al. 


(2010) concluded that there are simply not 
enough persistent systems for this to be the sole source 
of the gap: Even if every X-ray binary system that could 
contain a black hole were to contain a low-mass black 
hole, the gap would still not be fully populated. Nev¬ 
ertheless, it remains possible that a few of the missing 
low-mass black holes are lurking among transient X-ray 
sources. 

In this paper we present new optical photometry of 
V1408 Aql. We derive an improved orbital ephemeris for 
the system and model the mean orbital light curve us¬ 
ing our XRBinary light curve synthesis code, from which 
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Fig. 1.— Twenty nine nights of photometry of V1408 Aql. On 
the scale of this figure the individual nights are not resolved, only 
entire observing runs. The data in the four clumps on the l eft side 
of the figure were obtained by |Bayless et al. | poTT) and |Mason| 
|et al.| (J2012| , while the new data comprise the two clumps on the 
right side of the figure. The detached clump of bright points on 
the right side of the figure come from a single night, HJD 2456122 

we constrain the orbital inclination and the mass of the 
compact star. We find that the most likely mass for the 
compact star places it inside the gap in the mass distri¬ 
bution, although the range of possible masses is large. 
In section 2 we describe the observations and summa¬ 
rize the behavior of the light curve, and in section 3 we 
derive an updated orbital ephemeris. In section 4 we dis¬ 
cuss the models of the optical light curve and in section 
5 we discuss constraints on the distance to the source. 

TABLE 1 

Journal of Observations 


UT Date 

UTC Start 

Duration (hr) 

14-July-2012 

07:06 

3.6 

15-July-2012 

06:25 

4.5 

16-July-2012 

06:31 

4.2 

17-July-2012 

07:39 

3.1 

18-July-2012 

06:46 

4.0 

11-August-2012 

02:41 

4.7 

12-August-2012 

02:29 

4.5 

13-August-2012 

02:30 

4.0 

15-August-2012 

02:51 

4.9 


“Time resolution of all data is 10s. 


2. THE OPTICAL ORBITAL LIGHT CURVE 

We obtained new high-speed optical photometry of 
V1408 Aql with the Argos CCD photometer on t he 2.1m 
Otto Struve telesc ope at McDonald Observatory (|Nather| 
& Mukadam 20041. The photometer produces a sequence 
of consecutive CCD images, all with 10-second exposure 
times for the data reported here. We observed V1408 
Aql on five nights in July 2012 and four nights in August 
2012 for 3 - 5 hours per night (see Table fl|) . All the ob¬ 
servations were made through a broad BvR filter with 
a roughly-square passband from 4130 A to 7385 A. We 
reduced the data using standard IRAF routines and ex¬ 
tracted the brightness of V1408 Aql relative to the same 
two comparison stars used by Bayless et al. (2011). 

We combined the new photometry with the previously 
published photometry obtained using th e same cquip- 
rnent and reduced in a similar manner by|Bayless et ah| 
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(a) July 2012 lightcurves 
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(b) August 2012 lightcurves 

Fig. 2.— The light curves of V1408 Aql on five nights in July 
2012 (top) and four nights in August 2012 (bottom). 


(2011) and|Mason et al.| (2012). Altogether we now have 
29 nights of data extending from 2008 to 2013. The six- 
year light curve is plotted in Figure |T) On the scale of 
this figure the individual nights are not resolved, only 
entire observing runs. The data in the four clum ps on 
the left s ide o f the figure were ob tained by |Bayless et] 
al. (2011) and Mason et al. (2012), while the new data 
comprise the two clumps on the right side of the figure. 
The mean brightness of V1408 Aql varies by a factor of 
two. The night to night variations are of the order of one 
day, while the rapid variations or flickering are occasion¬ 
ally observed approximately every 1.1 ±0.2 hours. These 
long term variations might be caused by changes in the 
accretion rate of the secondary, which in turn affect the 
emission from the accretion disk. The detached clump of 
bright points on the right side of Figure [I] comes from a 
single night, 14 July 2012, when V1408 Aql was ~ 50% 
brighter than the immediately following night. Figure [2] 
shows just the nine new light curves. The variation in 
mean brightness from night to night is readily apparent 
as is the sinusoidal orbital modulation. Other sources 
similar to V1408 Aql show flares and bursts. For exam¬ 
ple, the low mass black hole LMXB GRO J0422±32 was 
found to have episodic gamma-ray and neutrino emission 
observed in the form of a hard power-law in the X-ray 


spectra (Vieyro et al. 


2012 ). 


The amplitude of the sinusoidal mo dulat ion is corre¬ 
lated with the mean brightness. Figure |3 (a) | shows a light 
curve from a nigh t when V1408 Aql was unusually bright 
and Figure |3(b) shows a light curve when it was unusu¬ 
ally faint. The sinusoidal modulation is clearly visible on 
both nights but its peak-to-peak amplitude was 29% on 
the bright night and only 22% on the faint night. Fig¬ 
ure [4] shows the amplitude of the sinusoidal modulation 
plotted against mean brightness from the six nights for 
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Fig. 3. — Two light curves of V1408 Aql. On June 06, 2008 (top) 
V1408 had a high mean brightness and the peak-to-peak amplitude 
of its orbital modulation was 29%. On June 25, 2009 (bottom) its 
mean brightness was lower and the amplitude of the modulation 
was 22%. 
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Fig. 4.— The peak-to-peak amplitude of the orbital modulation 
as a function of mean brightness from the six nights for which we 
have enough data to measure both accurately. The amplitude of 
the modulation increases as the mean brightness increases. 

which we have enough data to measure both accurately. 
There is a nearly linear correlation between the two with 
a slope dA/dB = 0.49, where A is the peak-to-peak am¬ 
plitude of the sinusoidal modulation and B is the mean 
brightness, both A and B are in units of relative inten¬ 
sity. 
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3. THE ORBITAL EPHEMERIS 

The night to night variations of the mean brightness 
tend to obscure the orbital modulation of the light curve, 
adding noise to measurements of the orbital period. The 
long term variations are on the order of days, or over 
twice as long as the orbital period. To mitigate this 
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Fig. 5.— The PDM periodogram for all our photometry of V1408 
Aql. The lowest minimum is the best orbital period. The points 
above the periodogram mark the previously publish ed or bital pe¬ 
riods from iT horstensenl (|1987|) , |Mason et al.| (|2012|l , and |Bayless| 
|et al.| ^Ollfr 

problem we multiplied the relative brightness on each 
night by a normalization factor to scale all the light 
curves to the same mean brightness, which we chose to 
be near the minimum mean brightness. There are a few 
caveats to mention. This works to first order because 
the long term variations are longer than the orbital pe¬ 
riod. And even though each light curve does not cover 
exactly the same phase, the normalization method for 
scaling should not be a problem. Since the normalization 
value used was less than 10% for more than 50% of the 
light curves used. We then measured the orbital period 
using the P hase Dispersion Minimization (PDM) peri¬ 
odogram (Stcllingwerf 1978|). Figure [5| shows the PDM 
periodogram along with the previously published periods 
and their error bars. The minimum of the periodogram 
is at a period of 0.388893(3) days, which is consistent 
with, but more accurate, than previously published peri¬ 


ods (Thorstensen||1987| Bayless et al.||2011 Mason et al 


20121. The improved orbital ephemeris is: 


T = HJD 2454621.829(4) + 0.388893(3 )E, (1) 

where T is the time of maximum flux and E is the orbit 
number. The five-year span of our data is too short for a 
meaningful constraint on the rate of change of the orbital 
period. 

Figure [6] shows the scaled data from all 29 nights folded 
at the orbital period. The sinusoidal orbital modulation 
is apparent, and so too is a large scatter about the mean 
orbital variation. The scatter is caused by rapid varia¬ 
tions - flickering - in the optical flux, which we presume is 
caused by rapid variations in t he flux from the accretion 
disk. Figure 8 of Mason et al. (2012) plots nightly light 
curves of V1408 Aql at a scale that displays the flicker¬ 
ing particularly well. The average orbital light curve is 
shown in Figure [6] 


4. ANALYSIS OF THE ORBITAL LIGHT CURVE 

We analyzed the orbital light curve of V1408 Aql with 
our XRbinary light curve synthesis program. All the 
models consisted of a black hole primary star surrounded 
by an accretion disk, plus a secondary star that fills 
its Roche lobe and is irradiated by flux from the disk. 
For our purposes the black hole is just a point source 
of gravity that emits no flux. The accretion disk is 
a cylindrically-symmetric, geometrically thin, optically 
thick, steady-state alpha-model disk. The outer radius 



Fig. 6.— Scaled photometry of V1408 Aql from all 29 nights 
folded at the orbital period. Each color represents a separate night. 
Phase zero is defined to be the maximum of the mean orbital light 
curve. The solid black line is the best-fit sine curve, not a model 
of the modulation. 

of the disk was set to the tidal truncation radius, which 
we take to be 0.9 times the mean radius of the primary 
star’s Roche lobe (Frank, King, & Raine 2002). We as¬ 
sume that the disk radiates like a black body. 

The flux emitted from a surface element of the sec¬ 
ondary star is prescribed to be F em i t = Fq + aFi rr , 
where F 0 is the flux the secondary star would emit in the 
absence of irradiation, Fi rr is the irradiating flux from 
the accretion disk, and a, which we call the “albedo,” 
is the fraction of Fi rr that is re-radiated instead of ab¬ 
sorbed into the structure of the secondary star. The 
intrinsic flux is calculated from the gravity darkening 
law F(j oc |.g| * * * 4 ^, where the local gravity is determined 
from the Roche geometry, and j3 is the temperature- 
dependent gravity-darkening coefficient, which we take 
from Claret (2000). We assume the emitted flux is fully 
thermalized, so the local effective temperature is given 
by aT* ff = F em it . Armed with a local effective 
gravity and local effective temperature, we adopt Ku- 
ruz solar-composition spectra for the local emitted flux 
if T e ff < 8000 K and blackbody spectra otherwise. The 
model does not include spots or other features on the 
surface of the secondary star. 

The free parameters of this model are 

• the masses of the primary and secondary stars. 

• the orbital inclination. 

• the inner radius and total luminosity of the accre¬ 
tion disk. 

• the intrinsic effective temperature of the (unirradi¬ 
ated) secondary star. 

• the albedo of the secondary star. 

There are also some additional parameters such as the 
time of phase zero and the resolution of the time steps 
that have no effect on the structure of the system. 

Perhaps surprisingly, our results do not depend 
strongly on the albedo of the secondary star. Changes in 
the value of the albedo are absorbed almost entirely into 
offsetting changes in the luminosity of the accretion disk. 
We have, therefore, fixed the albedo at 0.5. To estimate 
the inner radius of the disk we assumed the black hole 
has a spin of a* = 0.9 and used the sp in-ISCO relation 
in Figure 2 of ( McClintock et ak|2011 ). In fact, though, 
our results depend only weakly on the inner radius of 
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a M 2 = 0.8M© 
• M 2 = 1.0M, 
■ M 2 = 1.4M ( 
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Fig. 7.- Plot of the relative y 2 for the best-fitting synthetic light 
curves as a function of M\. The three sets of points correspond 
to the three masses of M 2 . We see that the best models are for a 
primary of 3Mq for any secondary mass. 


the disk. The inner radius does affect the fraction of 
disk flux emitted at optical wavelengths, but even here 
it is the outer disk radius that dominates the amount 
of optical flux. A numerical test showed that increasing 
the inner disk radius by a factor of three increased the 
preferred orbital inclination by only two degrees. The in¬ 
trinsic temperature of the secondary star is nearly irrele¬ 
vant since the optical flux is dominated by the accretion 
disk and the secondary’s heated face. Despite the large 
number of parameters needed to model the light curve, 
in the end only four parameters really matter: the mass 
of the primary star, the mass of the secondary star, the 
orbital inclination, and the disk luminosity. 

The light curve synthesis code assumes that the heat¬ 
ing of the secondary star is due to incident flux on it’s 
surface and that this flux is fully thermalized in the at¬ 
mosphere of the secondary star. Since the incident flux 
is fully thermalized, the details of the actual spectrum 
become irrelevant in favor of the total amount of inci¬ 
dent energy from the accretion disk. Therefore, the total 
flux depends on the total luminosity of the disk, not its 
spectral energy distribution. The disk luminosity is set 
via a separate input parameter in the code. 

The code assumes that the the spectrum emitted by 
an area element on the surface of the secondary is the 
same as the spectrum that would be emitted by an iso¬ 
lated star with the same effective temperature. The code 
does not account for the fact that the regions of the sec¬ 
ondary’s atmosphere which are heated by radiation can 
have a temperature inversion, which in turn can alter 
the absorption -line spectrum and produce chromospheri c 
emission lines (Barman et al. 2004 Wawrzyn et al. 2009). 
But since we are doing this analysis with optical photom- 
etry instead of spectroscopy, the form in which the code 
handles the incident radiation is not expected be the lim¬ 
iting factor in our results. 

We can eliminate inclinations greater than 65°. At in¬ 
clinations greater than i « 65° the accretion disk eclipses 
the secondary star. This is the case for all of the sec¬ 
ondary masses modeled (M 2 = 0.8 — 1.4M 0 ). Fig¬ 
ure 


8(a) 


shows an example of the distinctive eclipse pro¬ 
file produced when the accretion disk passes in front of 
the secondary star. No eclipse of this kind has been ob- 


servedF] To avoid an eclipse at 75° either the outer ra¬ 
dius ofthe accretion disk would have to be unrealistically 
small, about 50% the radius of the black hole’s Roche 
Lobe, or the mass ratio would have to be so low that the 
fits to the orbital light curve become unacceptably poor. 

Because so few parameters dominate the model, we 
were able to calculate a grid of models that covered the 
likely parameter space. The grid points in the primary 
mass were set at intervals 1 or 2 solar masses between 
2 Mq to 16 Mq (see the first column of Tables 2, 3, and 
4). The upper limit to the mass of a main-sequence sec- 
ondary star that just fills the Roche lobe is ~ 1.4 Mq 
(P atterson et al. 2005). We therefore chose grid points 
in secondary mass at 1.4M 0 , 1.0 M 0 , and 0.8 M 0 , the 
latter two corresponding to evolved secondaries. The grid 
points in inclination were typically at half degree inter¬ 
vals and in disk luminosity at intervals of roughly 10% 
of the luminosity. For the models with a secondary star 
mass of 1.4 M 0 we adopted a temperature of 6650 K, 
which is appropriate for a main-sequence star of that 
mass. For the models with secondary masses of 0.8 M 0 
and 1.0 M 0 we adopted temperatures of 4700 K and 
4900 K respectively, which are more appropriate for an 
evolved secondary. Again, though, the intrinsic temper¬ 
ature of the secondary has little effect on the synthetic 
light curve. The output at each grid point is a synthetic 
light curve and a value of relative % 2 for the fit of the syn¬ 
thetic light curve to the observed data. We report a rela¬ 
tive x 2 because the noise in the light curves is dominated 
by flickering noise, which is both highly correlated and 
variable in time, making measurement of the absolute x' 2 
intractable. The relative values of \ 2 are, though, ade¬ 
quate for determining the best-fit parameters and their 
standard deviations. Plots like the one shown in Figure 
[9] were used to pinpoint the inclination and disk lumi¬ 
nosity that yield the smallest relative x 2 for each pair 
of masses. In that Figure, the red zones have the high¬ 
est values of x 2 , while the blue regions have the lowest 
values and correspond to better fits to the data. Once 
the masses have been chosen, the inclination is highly 
constrained, typically to within 0.5°, the disk luminosity 
less so, but still typically to within 30%. 

The results are given in Tables [2j [3j and |4j which list 
the relative x 2 for the the best-fit values ofthe inclina¬ 
tion and disk luminosity for each pair of primary and 
secondary masses. Figure [7] intercompares the models 
by plotting the best-fit values of x 2 against the mass 
of the primary star. The lowest values of x 2 occur for 
M 2 = 0.8 and 1.0 Mq, and for Mi between 2 and 5M 0 . 
The 90% upper bound on the mass of the black hole is 
6.2 M 0 . We are unable to place a useful lower limit on 
the mass of the black hole, but masses in the range of the 
known neutron star masses are not excluded. The very 
best fit occurs for Mi = 3.0 Mq, M 2 = 1.0 M 0 , and an 
orbital inclination i = 12.75°. The preferred mass of the 
secondary in our models is perfect ly consistent with the 
results found by Hakala et al. (2014). Our models show a 
90% upper bound of 6.2 Mq, and "a 90% lower bound that 
is outside the reasonable mass for a black hole primary. 

4 If the disk eclipses the secondary star, the secondary star also 
eclipses the disk 1/2 orbit later, but a shallow disk eclipse can be 
difficult to discern because it is superimposed on the minimum of 
the sinusoidal orbital variation. 
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(a) The best-fitting synthetic light curve for an 
inclination of 75 °. The dip near phase zero is 
caused by an eclipse of the secondary star by the 
accretion disk. 



Phase 


(b) The very best fitting synthetic light curve. The 
parameters of the model are Mi = 3.0 Mq, M 2 = 
1.0M o , and i = 12.75° 


Fig. 8. — Two best-fit synthetic light curves overplotted on the 
mean orbital light curve of V1408 Aql. There is a difference be¬ 
tween the average light curve and the model fits at phase = 1.0. 
This is likely due to insufficient coverage at this phase. However, 
additional photometry is needed to confirm this result. 


The lcr deviation from the best fit of Mi = 3.0 Mq is of 
±2.5 Mq. Even though we can not constrain the lower 
mass based solely on our modeling, we adopt a lower 
bound of Mi = 2.0 Mq from the assumption that the 
primary is a black hole and not a neutron star. Black 
holes are expected to form only when the mass of the 
compact object is too high to become a neutron star. 
Some ways in which a black hole can be created include: 
the collapse of the core of a star during a core-collapse 
supernova, a merger of two neutron stars, and the col¬ 
lapse of a neutron star that is accreting mass. In all these 
scenarios the black hole will not form, unless the mass 
is greater than the highest possible mass for a neutron 
star. Currently, the observed upper mass limit of ncu- 
tron stars is 2.0 Mq ( Demorest et al.|[2010 1. Figure 8(b) 


shows the synthetic light curve tor this model overplotted 
on the mean optical light curve. The correlated errors in 
the light curves and imperfect model physics broaden the 
range of permitted black hole masses beyond the formal 
90% upper limit, possibly even to much higher masses. 
The best-fit masses are, however, unchanged by these 
considerations. 


5. DISCUSSION 

The fits to the orbital light curve of V1408 Aql favor 
low orbital inclinations, much lower than the inclinations 
usually adopted in analyses of its X-ray spectral energy 



5 6 7 8 9 10 11 12 13 14 

Disk Inclination (degrees) 


Fig. 9. — The contour plot represents various models with differ¬ 
ent inclinations and disk luminosities for the case of a 2 Mq primary 
and a 0.8 Mq secondary. The red regions represent the worst fits 
and the blue regions represent the lowest \ 2 an d therefore the best 
fits to the optical light curve. The red cross marks the spot of the 
best fit with an inclination of i = 11 ° and a luminosity of L = 
7.50X10 37 ergs/s. 


distribution. According to Maitra et al. (2013), for exam¬ 
ple, fits to the X-ray spectrum of 4U 1957 • 115 require 
an orbital inclination near 75° to obtain the canonical 
X-ray spectral hardening: factor hd — T co i or /T e ff = 1.7 
(Davis et al. 2005 20061 (this quantity is often called 
the color correction factor and denoted by f c ). Lower 
inclinations yield values for hd that are greater than 1.7. 
To rigidly limit the color correction factor to values near 
1.7 is, however, unwarranted. Authors such as|Nowak et| 


al. (2008) preferred high inclinations near 75°, which in 
combination with a large distance, low black hole mass 
and high accretion rate reconcile the high temperature 
and low normalization obtained from their models. For 
some of their models they fixed the inclination to 75° 
in order to be cons istent with th e interpretation of opti¬ 
cal observations of Hakala et al. (1999). There is ample 
theoretical and observational evidence that color correc¬ 
tion factors can be much larger than 1.7, Cygnus X-l 
shows o ne of t he hig hest color correction factors of up to 
/ c ~ 5 (Reynolds et al. 2013). Color correction factors 
are different for various X-ray sources, and beyond that, 
an individual source can have different color correction 


factors (Merloni et al. 

2000 

Dunn et al. 

2011 

Salvesen 

et al. 2013|T 


tral energy distribution, we note that color correction 
factors larger than 1.7 do permit reasonable system pa¬ 
rameters for low orbital inclinations. Although the typ¬ 
ical value for the_spectral hardening factor is f c = 1.7, 


Maitra et al. (2013) provided a val ue o f f c = 2.0 — 2.2 
for their preferred model. Nowak et al. (20121 fit the 
spectral energy distribution of of 4U I957±115 with sev¬ 
eral sets of models, one of which was the eqpair model 
with color correction factors raging from f c = 1.7 — 3.3. 
While the fits of the eqpair model yielded acceptable 
values for y; 2 , they required a low valu e for the n ormal- 

~' (12012 1 


Nowak et al. 


ization factor, N eqp = 1.926 x 10 -4 . 
interpreted the low normalization factor as evidence tor 
a small inner disk radius and, consequently, evidence for 
a rapidly spinning black hole. While the fits all assumed 
an orbital inclination of 75°, the p aramet ers of the fits 
are highly degenerate. Nowak et al. (2012) give a scaling 



























































TABLE 2 

Models with a Secondary Star Mass of O.8M 0 
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Mi 

(Mq) 

a 

(AU) 

Disk 

(a) 

i 

(deg) 

Ld 

(10 38 erg s -1 ) 

Lo/LEdd 

Relative x 2 
(xlO 5 ) 

2 

0.0145 

0.414 

11.0 

0.75 

0.30 

1.121 

3 

0.0160 

0.446 

13.5 

0.72 

0.19 

1.120 

4 

0.0173 

0.469 

16.0 

0.71 

0.14 

1.122 

5 

0.0184 

0.486 

18.5 

0.78 

0.12 

1.121 

6 

0.0194 

0.499 

21.0 

0.70 

0.09 

1.124 

7 

0.0204 

0.511 

23.5 

0.72 

0.08 

1.130 

8 

0.0212 

0.520 

26.0 

0.75 

0.07 

1.140 

10 

0.0227 

0.536 

31.0 

0.90 

0.07 

1.156 

12 

0.0240 

0.548 

36.0 

1.30 

0.09 

1.171 

14 

0.0252 

0.558 

41.0 

1.60 

0.09 

1.177 

16 

0.0263 

0.567 

46.0 

2.90 

0.14 

1.162 


Note: a is the separation of the stars, Roisk is the outer radius of the disk in units of a, Ljj is the luminosity of the disk, and L^dd is the 
Eddington luminosity. 


TABLE 3 

Models with a Secondary Star Mass of 1.0M 0 


Mi 

(Mq) 

a 

(AU) 

R-Disk 

(a) 

i 

(deg) 

Ld 

(10 38 erg s" 1 ) 

Lo/LEdd 

Relative \ 2 
(xlO 5 ) 

2 

0.0148 

0.396 

10.50 

1.4 

0.56 

1.121 

3 

0.0163 

0.429 

12.75 

1.3 

0.34 

1.118 

4 

0.0176 

0.451 

15.00 

1.4 

0.28 

1.120 

5 

0.0187 

0.469 

17.25 

1.4 

0.22 

1.121 

6 

0.0196 

0.482 

19.50 

1.6 

0.21 

1.123 

7 

0.0205 

0.494 

21.75 

1.6 

0.18 

1.125 

8 

0.0214 

0.504 

24.00 

1.6 

0.16 

1.130 

10 

0.0228 

0.520 

28.50 

1.9 

0.15 

1.150 

12 

0.0241 

0.533 

32.50 

1.7 

0.11 

1.165 

14 

0.0253 

0.544 

36.50 

1.7 

0.10 

1.177 

16 

0.0264 

0.552 

41.00 

2.8 

0.14 

1.182 


See note to Table [2] 

relation for other system parameters: 


N = 

1 y eqp 


( M 1 
\1M q 



( 2 ) 


where D is the distance and f c is the color correction 
factor. With M\ = 3M 0 , i = 12.75°, and N eqp = 1.926x 
10 -4 , this becomes 


Df c = 213.5. 


(3) 


If one insists on f c = 1.7, the distance would be D = 
74 kpc, placing V1408 Aql uncomfortably far outside the 
Galaxy. If, instead, we set f c = 2.3, the distance becomes 
40 kpc; and if f c = 3.3, it drops to 20 kpc, either of which 
would place V1408 Aql in the halo of the Galaxy. 

This range of distances agrees with the results of[R us- 
sell et al. (2010). If V1408 Aql contains a black hole, it 
would need to De at a distance between 22 and 40 kpc 
to lie in the same region of the (Lx, Lqpt) diagram as 
other black hole X-ray binaries in their soft states, with 
the l a rger dis t ance s omewhat preferred (see Figure 8 in 
Russell et al. (2010)). At the larger distance the X-ray 
luminosity at 2-10 keV would be ~ 10 38 erg s _1 . This 
agrees with the large, distance-independent disk lumi¬ 
nosity that is required to heat the secondary star in our 
models of the orbital light curve. 


Finally, we have found a correlation between the mean 
brightness of V1408 Aql and the amplitude of the sinu¬ 
soidal orbital modulation: The amplitude inc reases as 
the mean brightness increases (see Figure [4]). Russell et 


al. (2010) found a correlation between the mean optical 


and mean X-ray fluxes from V1408 Aql: The mean X-ray 
flux increases as the mean optical flux increases. Taken 
together, these two correlations imply that the amplitude 
of the sinusoidal orbital modulation increases as the X- 
ray flux increases. This is the correlation expected from 
a model in which the orbital modulation is caused by 
the heated face of the secondary star. The X-ray flux is 
a good proxy for the total disk luminosity. As the X- 
ray flux increases, the disk becomes more luminous and 
heats the face of the secondary star towards the disk to 
a higher temperature, increasing the amplitude of the 
orbital variation. 


6. SUMMARY 

We have presented new optical high-speed photometry 
of V1408 Aql from nine nights in 2012 July and August. 
The optical light curve continues to display a nearly- 
sinusoidal orbital modulation along with night-to-night 
variations of the mean brightness. We combined the new 
photometry with our previously-published photometry to 
derive a more accurate orbital period and mean orbital 































TABLE 4 

Models with a Secondary Star Mass of 1.4M® 


Mi 

(Mq) 

a 

(AU) 

Disk 

(a) 

i 

(deg) 

Ld 

(10 38 erg s -1 ) 

Lo/LEdd 

Relative x 2 
(xlO 5 ) 

2 

0.0154 

0.369 

10.0 

1.2 

0.48 

1.129 

3 

0.0168 

0.401 

12.0 

1.7 

0.45 

1.127 

4 

0.0180 

0.425 

14.0 

1.9 

0.38 

1.127 

5 

0.0191 

0.442 

16.0 

2.3 

0.37 

1.129 

6 

0.0200 

0.457 

18.0 

2.8 

0.37 

1.134 

7 

0.0209 

0.469 

20.0 

3.0 

0.34 

1.140 

8 

0.0217 

0.479 

22.0 

3.4 

0.34 

1.145 

10 

0.0231 

0.496 

26.0 

4.2 

0.33 

1.161 

12 

0.0244 

0.509 

30.0 

5.0 

0.33 

1.183 

14 

0.0256 

0.520 

34.0 

6.3 

0.36 

1.208 

16 

0.0266 

0.530 

38.0 

8.9 

0.44 

1.224 


See note to Table 01 

light curve, and to better define the night-to-night vari¬ 
ations. We find that the amplitude of the orbital mod¬ 
ulation is strongly correlated with the mean brightness, 
dA/dB = 0.49, where A is the peak-to-peak amplitude 
of the sinusoidal modulation and B is the mean bright¬ 
ness, both A and B are in units of relative intensity. The 
relative amplitude of the orbital modulation rises from 
23% when V1408 Aql is at the minimum of the observed 
range of its brightness to 29% at the maximum of the 
range. We attribute the changes in mean brightness to 
changes in the luminosity of the accretion disk around 
the black hole. 

After scaling all the nightly light curves to the same 
mean brightness, we derived a more accurate orbital pe¬ 
riod, 0.388893(3) days, and mean orbital light curve, 


shown on Figure 8(b) The mean orbital light curve is 


consistent with a model in which the orbital modulation 
is caused entirely by the changing aspect of the heated 
face of the secondary star as it revolves around the black 
hole. Fits of synthetic orbital light curves based on this 
model to the observed light curve favor low orbital incli¬ 
nations and low black hole masses, the best fit occuring 
for Mi = 3.0 Mq, M 2 = 1.0 M 0 , and i = 12.75°. The up¬ 
per bound to the mass of the black hole is 6.2 Mq with a 
90% probability, although uncertainties in the data and 
the models allow higher masses, possibly much higher 
masses. Orbital inclinations higher than about 65° are 
strongly disfavored by the lack of eclipses. 

The low orbital inclinations we have found are com¬ 
patible with previous analyses of the X-ray spectral dis¬ 
tribution of V1408 Aql if the color correction factor is 
somewhat larger that the value typically adopted for the 
analyses. If the distance to V1408 Aql is 40 kpc, the 
color correction factor must be increased to 2.3; and if 
the distance is 20 kpc, the color correction factor is 3.3. 

In conclusion, the compact star in V1408 Aql a viable 
candidate for a black hole whose mass lies within the gap 
in the distribution of compact star masses. 
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